The Hydraulic Jump In Two 
Dimensions 



Abhishek Mukherjee , Amaresh Datta 
Jayanta K Bhattacharjee 2 

Department of Physics, Indian Institute of Technology Bombay-400476 
2 Department of Theoretical Physics, 
Indian Association for the Cultivation of Science, 
Jadavpur, Calcutta - 700032, India, 
e - mail: tpjkb@iacs.ernet.in 



PACS : 47.10. +g ; 47.20 Gv 

Abstract 

The first quantitative calculation of the position of the hydraulic 
jump was carried out a few years ago by Bohr et.al. Since this is 
the only calculation of the jump, we have analysed it from a slightly 
different point of view. Our results are similar to those of Bohr et.al. 
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If a vertical jet of fluid impinges on a horizontal surface and spreads out 
radially, then at a certain radius one observes a sudden jump in the height 
of the fluid[l-4].It is a familiar observation, seen in the kitchen sink everyday, 
but it has been a difficult system to deal with. The question that is normal 
to ask is that given the volumetric flow rate (or the Reynolds number of 
the impinging flow) can we predict the position and magnitude of the jump 
which is known as the hydraulic jump. The first attempt to provide a concrete 
answer to this question was undertaken by Bohr,Dimon and Putkaradge[5] 
and took place only a few years ago. 

Continuity of mass flux and momentum flux is possible at any radius and 
thus one more condition is needed to determine the radius of the jump. This 
could be the energy flux but that is generally difficult to handle. It was shown 
by Bohr et.al. that introducing viscosity does help in setting up a criterion 
for finding the position of the jump. Since that was the first attempt to set 
up a quantitative calculation, we believe that it is worthwhile to look at the 
problem from a simpler standpoint. 

The central simplification in the calculation, a simplification also used by 
Bohr et.al., is the assumption that the radial velocity Vis greater than the 
vertical velocity 'w', while the variations with respect to the vertical direction 
(ie.V) is far more important than the radial variations .With this in mind, 
the equation of continuity becomes: 

+ Q- z ( Wr ) = 

du u dw ,^ 
dr r dz 

The momentum flow equuation(Navier Stokes equation), for the radial com- 
ponents: 

du du 1 dp d 2 u 
dr dz p dr dz 2 
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If Q is the volumetric flow rate, 

rh(r) 

Q = 2nr / u(r,z)dz, (3) 
Jo 

where h(r) is the height of the fluid layer at distance r.The boundary con- 
dition at z=0 is that the velocity variations is constant with the plate 
ie.u=0=w.On the surface ie.z=h(r)the stress is zero ie. |^ = O.We now 
make the reasonable Ausatz that 

u(r,z) = U(r)f(r,) (4) 

where r\ is the scaled variable 77 = j4^r. The boundary conditions decree that 
f(0)=0 ;f(l)=l and /(/)(1) = O.From Eqn(3) we have 

Q = 2nrU(r)h(r) [ f(rj)dr) 
Jo 

= 27vrU(r)h(r)C (5) 
where C = Jq 1 f {r])dr}. From Eqn(l) we have 

u(r z) = _9_M_ (6) 

U[r,Z) 2nCrh(r) W 

and 

w(r,z) = U(r)h'(r)rif(r,) (7) 
Turning to Eqn(2), we see 

or] or 
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and 



leading to 



dz h dr) 
1 d 2 f 

- "W < 10 » 

p or dr 

W'f- - + f (12) 

At this point, we need to make a statement about the vertical profile f(rj).A 
common situation (constant pressure gradient)is one where the profile is 
parabolic. This implies (with the boundary conditions f(0)=0,f(l)=l and 
/(/)(1) = 0)that f(rj) = 2(77) - (??) 2 .this yields /(/)(/) = -2 and we get 
on the surface (ie.77 = 1) 

UU' = -gti - ^ (13) 

which is the main equation for considering the phenomenon of hydraulic 
jump. Bohr et.al. arrive at a similar equation. using the continuity condi- 
tions as expressed by Eqn(5), U(r)h(r) = J^-, we can write, U'{r)h{r) + 
U{r)h'{r) = —-^2 and this helps us to write Eqn(13) 

3Q U(r) h'U 2 2vU 

-gh 



47rr 2 h(r) h " h 2 

dh, U 2 , , 3Q 2uM 
-iz(9--r) = [ 



dr h Anr 2 h h 

h— = \j\Mxi hi (14) 
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the central equation of our paper. We can write this as a set of two equqtions 

ldh 2 3Q 2v 3Q 

2 dr Anr 2 h Aurh 

dr U 2 ,3Q 2 1 

Tr = g -T = 9 - { A^V (15) 

This system has a fixed point at r=R,U = U ,h = h with 

U% = gh (R) 
M«) = |f (16) 

Writing h = %H,r = ^^r,u = J~(®?)U. U 2 = H and H = i? 2 .Thus R 

" V (■?■) 1 " 3 

at the fixed point is: R = ^(f) 2 (^-)^Q~' § (v)«g~«. THus in unsealed coor- 

5 3 1 

dinates, r^atfixedpoint) = (0.99)gs (z/)~ig~i .It is straightforward to check 
that the fixed point is astable spiral. 

We now return to Eqn(14) and examine the limits r ^> R and r <C i?.For 
r i?,Eqn(14) simplifies to 

dh ^ 3Qv 
dr 2irgrh 3 

which leads to 

6g, ln ro ]1/4 (lg) 
n g r 

for r ^> R but r <C r .the outer cutoff is r beyond which this formula is 
not valid. We see that h(r) has a very weak r-dependence in this range. for 
r <C i?,on the other hand, we have 

Hh 3Q 3Q 

u/t _ , 4vrr/i 47rr 2 



which integrates to 



rh = constant (20) 
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In the inner region the height drops off as h(r) ~ - and in the outer region it 
is approximartely constant. The dividing point of the flow is at r=R,which is 
where the jump takes place so that falling h ~ \ profile can meet the nearly 
constant h ~ profile. 

This leads to the jump radius being 

R= ^y/s^^s/s -1/8 (21) 

This relation is very similar to that found by Bohr et.al. The difference in the 
prefactor comes from the difference between the choice of the z-dependence 
of the radial velocity.Bohr et.al. chose to work with a z-averaged velocity 
while we have demonstrated the jump by studying the radial flow velocity at 
the free surface. 
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